EbN0_db_nc = 0:11;
Ber_nc = [0.0786530    0.0562721    0.0375477    0.0228728    0.0124905    0.0059547    0.0023887    0.0007745  ..
0.0001925  0.0000326  0.0000037    0.0000002];

EbN0_db_74 = 0:7;
Ber_74 = [0.0441516    0.0241727    0.0113311    0.0044057    0.0013594    0.0003139    0.0000513    0.0000046];

EbN0_db_157 = 0:6;
Ber_157 = [0.0317736    0.0138362    0.0047395    0.0012093    0.0002142    0.0000252    0.0000014];

EbN0_db_conv = 0:5;
Ber_conv = [0.020632    0.0078374    0.002305    0.000494    0.0000816    0.0000104];

plot(EbN0_db_nc, log10(Ber_nc), EbN0_db_74, log10(Ber_74), ..
  EbN0_db_157, log10(Ber_157), EbN0_db_conv, log10(Ber_conv))
//plot(EbN0_db, log10(0.5*erfc(sqrt(EbN0))))

